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Abstract 



It is shown that a normalisable probability density can be defined for the entire 
complex plane in the modified de Broglie-Bohm quantum mechanics, which gives 
complex quantum trajectories. This work is in continuation of a previous one that 
defined a conserved probability for most of the regions in the complex space in terms 
of a trajectory integral, indicating a dynamical origin of quantum probability. There 
it was also shown that the quantum trajectories obtained are the same characteristic 
curves that propagate information about the conserved probability density. Though 
the probability density we now adopt for those regions left out in the previous work is 
not conserved locally, the net source of probability for such regions is seen to be zero 
in the example considered, allowing to make the total probability conserved. The 
new combined probability density agrees with the Born's probability everywhere on 
the real line, as required. A major fall out of the present scheme is that it explains 
why in the classical limit the imaginary parts of trajectories are not observed even 
indirectly and particles are confined close to the real line. 
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1 Introduction 



Complex quantum trajectories were first obtained pQ by modifying the de 
Broglie's guiding wave approach to quantum mechanics. Here, trajectories 
were drawn for the cases of harmonic oscillator, potential step, wave packets 
etc. For getting this trajectory representation, first we substitute \l/ = e %s / h in 
the Schrodinger equation, which gives the quantum Hamilton-Jacobi equation 
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Then we postulate an equation of motion for the particle, similar to that used 
by de Broglie: 

. = os_ h i_a^ 

dx i $ dx 



The quantum trajectories x(t) were found by integrating this equation with 
respect to time; in general, they lie in a complex x-plane, with x = x r + ixt. 
This results in a modified version of the de Broglie-Bohm (dBB) quantum 
mechanics [415116117] . Eq. (J5J) was used by Leacock and Padgett [8 J to obtain 
eigenvalues in many bound state problems, without having to solve the corre- 
sponding Schrodinger equation. 

It shall be noted that the canonical momentum is not always the mechani- 
cal momentum, even in classical mechanics with Cartesian coordinates. For 
instance, when there are velocity-dependent potentials, the canonical and me- 
chanical momenta are different. In the case of a charged particle in a magnetic 
field with the vector potential A , the mechanical momentum may be written 
as 

mx = V5- eA/c. (3) 



This equation in Cartesian coordinates shall thus be the equation of motion 
we adopt for charged particles in higher dimensions with electromagnetic field. 

The Floyd-Faraggi-Matone (FFM) trajectory representation [9|10|I11|I12|[T3"] is 
another modified dBB version but with real trajectories and is based on a 
generalised Hamilton- Jacobi equation equivalent to that used in dBB. But 
this representation differs from dBB quantum mechanics mainly in the use of 
the equation of motion. Here, for stationarity the equation of motion for the 
trajectory time t, relative to its constant coordinate r, is given as a function 
of x by 

t- T = OW/dE (4) 



where W is the Hamilton's characteristic function appearing in the generalised 
Hamilton- Jacobi equation and E is the energy. Carroll [11] finds that for sta- 
tionarity the above Jacobi's theorem is valid, for W is a Legendre transform 



of Hamilton's principal function. Floyd [12] notes that as Jacobi's theorem 
also determines the equation of motion in classical mechanics, it is universal 
transcending across the division between classical and quantum mechanics. In 
this way, FFM claims to be a deterministic theory. 

The complex function S, which may be called the complex Hamilton's princi- 
pal function in the present modified de Broglian mechanics, and the complex 
QHJE (TjQ) itself are quite different from the corresponding entities in dBB 
or FFM representations. Moreover, probability is an integral part of the new 
formulation and one can also consider equation (T5]) as propagating probabil- 
ity densities. In fact, it was shown in a previous work [T5] that the quan- 
tum trajectories obtained in this scheme are the same as the characteristic 
curves propagating information about a conserved probability density. Hence 
we continue to use ([2]), with appropriate modifications as mentioned above 
when necessary, as the equation of motion in this modified dBB formulation. 

The complex trajectory approach gives the paths in the n = 1 harmonic 
oscillator as shown in Fig. 1. In this case, the paths are given by (a 2 x 2 — 
a 2 x 2 — l) 2 + 4a 4 x 2 x 2 = A 2 , a constant for each path. Here A is real and 
positive, and a 2 has the usual definition, equal to muo/h. It is interesting to 
note that these figures are the famous Cassinian ovals. The Jacobi lemniscate, 
which is the special case of these ovals, corresponds to A = 1 . 




Fig. 1. The complex trajectories in the n = 1 harmonic oscillator case, drawn for 
various values of A and with axes X r = ax r and Xi = axi. These are the famous 
Cassinian ovals. The lemniscate corresponds to A = 1. 

Recently, the prospects of a dynamical explanation of quantum probability 
in this scenario was explored [H]. We have shown that in this modified de 
Broglian mechanics, the Born's probability along the real axis can be obtained 
as an exponential function of a line integral along the real axis. Similarly in 
the extended x r Xj-plane, a probability density is obtained as the exponential 



of a trajectory integral. In both these cases, integrands involve the particle's 
complex velocities, indicating a dynamical origin of quantum probability. In 
the extended case, we obtain a conserved probability, which agrees with the 
boundary condition (Born's rule) along the real axis in most regions. However, 
for other regions (e.g., the subnests inside the lemniscate in the n — 1 harmonic 
oscillator case, which do not enclose any poles of x), such a definition could not 
be obtained. In fact, a probability density was not proposed for such regions 
in [H] and it was considered a missing link in the formalism. 



Another important question faced by this trajectory representation is why the 
complex extended motion is not appreciable in the classical limit. Our every- 
day experience is that particle trajectories are in real space, obeying classical 
rules. In this letter, we study the probability distribution of particles in the 
complex plane in detail and adopt for those regions left out in [13] an alter- 
native definition, first proposed in [TofTo] as the probability density for the 
entire plane. This density satisfies the boundary condition on the real axis, 
though it is not a conserved one. We show that the combined probability for 
the entire extended plane does not diverge and there is no difficulty in nor- 
malizing it, contrary to what happens in the above case. Nonconservation of 
probability points to the prospects of particle creation and destruction, a fun- 
damental feature of quantum phenomena. Most importantly, the new scheme 
reveals that large part of the probability for a classical harmonic oscillator lies 
extremely close to the real line, thereby explaining why the complex extended 
motion does not leave any imprint in the classical limit. 

We restrict ourselves to one particle stationary states in 1-dimension for sim- 
plicity. The paper is organized as follows. In the next section, we review the 
problem with probability density faced in |14j . The third section presents the 
alternative definition of probability density adopted for the subnests and ob- 
tains a form of continuity equation obeyed by it. In section 4, the classical 
limit of a harmonic oscillator in the light of our combined probability density 
and its implication for complex trajectories are discussed. Section 5 comprises 
our conclusion. 



2 Probability from velocity field 



As in the case of dBB quantum mechanics, the present modified de Broglie- 
Bohm mechanics too is constructed in such a way that it agrees with the results 
of standard quantum mechanics |TJ. This is achieved by accepting the Born's 
probability axiom along the real line in these representations. However, we note 
that one of the challenges before such a quantum trajectory representation 
is to explain this probability axiom. In the standard dBB approach, there 
were several attempts to obtain the \I>*\I> probability distribution from more 



fundamental assumptions [6]. (The FFM trajectory representation does not 
involve probability and is considered a deterministic description.) 

In addition to that defined along the real line, it is desirable to have a proba- 
bility density defined everywhere in the extended complex x-plane. An earlier 
attempt made in [T7] to define such a density was to write p(x) = \l/(x)\l/(x), 
where ty(x) = ty*(x*), with x complex. With the help of time-dependent 
Schrodinger equation, the author shows that, in general, dp/dt ^ j'(x,t). 
This arguably leads to nonconservation of probability along trajectories. But 
it shall be noted that this negative result is based on the choices made in [T7] 
for the probability density and flux. Moreover, this definition leads to com- 
plex probability off the real axis, which is undesirable. Another proposal is to 
define probability as ty*(x)^(x) itself [T5iri6] . Though this has the advantage 
of being real everywhere, it is not conserved at any point in the extended 
plane. In addition, this diverges for large values of Xi in many cases and is not 
generally a normalisable probability. 

Using a different approach, it was shown in [2] that the complex trajectory 
representation is capable of explaining the quantum probability as originating 
from dynamics itself. Here, the Born's probability density to find the particle 
on the real axis around some point at x — x r o was obtained as 

r- / XidXr , (5) 



where the integral is taken along the real line. In addition, an extended prob- 
ability density p(x r ,Xi) in the x r Xj-plane for stationary states was proposed 
in [2] and showed, with the aid of complex-extended Schrodinger equation, 
that the continuity equation follows from it. The proposal was that if po, the 
extended probability density at some point (x r0 , Xjo) is given, then p(x r , Xi) at 
another point that lies on the trajectory which passes through (x r0 ,Xio) is 



p(x r (t),Xi(t)) =p exp 



-4 }, n 



h 

to 



I Im i-mx 2 + V{x)\dt' 



(6) 



Here, the integral is taken along the trajectory [x r (t'),Xi(t')] which passes 
through (x r o,Xio). The continuity equation was shown to follow from this 
axiom by using the extended version of the Schrodinger equation. While eval- 
uating p with the help of © above, one should specify p at (x r o,%io) an d if 
we choose this point as (x r0 , 0), the point of crossing of the trajectory on the 
real line, then p Q can assume the Born's value P(x r0 ). 

However, we may note that in some regions of the complex plane, there can be 
disagreement between the values of p and P(x r ) at the points of crossing of the 



trajectories on the real line. This happens where the complex trajectories do 
not enclose any poles of x, described as 'subnests' in jT]. Stated more clearly, 
the problem here is that as a particle trajectory is traversed in these regions, 
even while the probability p at one point of crossing x r0 on the real axis agrees 
with P(x r0 ), at the other point, say the point x r \ where the trajectory again 
crosses the real line, the probability calculated according to (Q happens to be 
different from that of P(x r i). 

On the other hand, if we directly solve the continuity equation for stationary 
states, it is easy to see that its solution and the probability density given 
by the trajectory integral ([H]) give identical results for regions outside the 
subnests [H]. But for the subnests, the boundary condition overdetermines 
the problem and we are unable to find a solution that agrees with the Born's 
rule everywhere on the real line; i.e., there is complete agreement between the 
extended probability density p and the Born's probability density only in the 
regions outside the subnests. 



3 Probability inside the subnests 



Given this situation, one can ask whether it is possible to find a probability 
density for the subnests that can agree with the Born's rule on the real line, 
even if it is not conserved in the region. A natural choice for such a definition 
is the extended probability p' = ty*(x)^(x), suggested in [T5][T6] . A trajectory 
integral form for p', similar to that in Eq. (EJ) above, was proposed in [18J. For 
stationary states, this can be written as 



p'(x r ,Xi) ocp exp 
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(7) 



We propose to adopt this as the probability density for such subnests. Here 
also, the integral is to be evaluated over the trajectory of the particle. The 
absence of the potential term V(x) in the integrand marks this definition from 
that in ([6]). Therefore this distribution will not be conserved. Instead of the 
continuity equation, we get 

dp' d(p'x r ) d(p'xi) 4 .. . T , Tr/ .-, 



Thus there are sources and sinks for probability in the subnests. However, as 
seen below, the net source of probability in this region for the n = 1 harmonic 
oscillator is zero, indicating that the total probability for the subnested region 
can remain conserved. 



Thus instead of being 'nonviable', the trajectories in this region reveal an im- 
portant feature of quantum phenomena, namely nonconservation of particles. 
This indicates creation and destruction of particles in such regions. Hence 
we can anticipate that the representation is capable of allowing a smooth 
transition to quantum field theory. It is interesting to note that the Floyd- 
Farraggi-Matone trajectory representation too allows creation and destruction 
of particles [13]. 



4 Probability and classical limit 



The extended probability density for the n — 1 harmonic oscillator in the 
region inside the lemniscate (for A < 1 and x r > 0), computed using the 
trajectory integral approach in ([7]) is shown overlapped with the 'leaf-shaped' 
surface \I/*\I/ in this case, in Fig. 2. 





Fig. 2. The extended probability density p(x r ,Xi) for the A < 1, x r > region 
of the n = 1 harmonic oscillator, evaluated along the trajectories, overlapped with 
the extended v]/*^ probability density. This probability density does not obey the 
continuity equation, but the total probability will be conserved. Also it agrees with 
the Born rule along the real line. 



Using the above definition p' = \|>*\|/ and the expression for x in the n 
case, we shall see that 



d(p'x r ) d(p'xi) 



dx T 



+ 



dxi 



OC e q2(:e2 ^Hx^f + ^Xj 



(9) 



The quantity on the right hand side is the density distribution of the 'prob- 
ability source'. It is positive for the first and third quadrants and negative 
for the second and fourth ones. When the entire region inside the lemniscate 
is considered, the net source of probability is seen to be zero. Thus the total 
extended probability can be normalized for the n — 1 harmonic oscillator. 

It is found that the probability in the region inside the lemniscate is sub- 
stantial; 43.25% of the total probability lies inside this region for the n — 1 
harmonic oscillator. The maximum value of Xi for the lemniscate (a measure of 
its width) in this case is x™~ ax = X^^/a = 0.4858/a. This explains how classi- 
cal particles are confined close to the real line. For instance, consider a classical 
harmonic oscillator of mass m ~ 1 kg and frequency ~ 1 Hz in the n = 1 state. 
The maximum value of x, for its lemniscate is x™ a;c ~ 10 ~ 17 / y/muo in units 
of metres. Thus 43.25% of the total probability in this case lies within the 
lemniscate of width ~ 10~ 17 m. It may also be noted that since the extended 
probability for A > 1 decreases fast, most of the probability outside also lies 
close to the real line. 

For the higher energy eigenstates of the harmonic oscillator, the lemniscates 
are seen to be narrower than that of the low energy ones. For example, in the 
n = 2 case, the complex paths are described by 



(x; + x;y - 5{x; - xt) + 



25 



i 2 



(X% + Xf) = constant. 



(10) 



These are shown in Fig. ([3]) . The X™ ax for the 3-fold lemniscate in this case is 
0.4125 and hence its x™ ax is less than that in the previous case. For n — 3, the 
complex paths are described by a more detailed expression, and are shown in 
Fig. (HI). The X™ ax for this 4-fold lemniscate is found numerically as ~ 0.39, 
which is again smaller than that in the previous cases. Thus it can be assured 
that for a classical oscillator of any energy and having mass 1 kg. , the width of 
the region, where the large part of probability lies, is less than or of the order of 
10~ 17 / \/uJq m. This indicates that the probable imaginary part of position for 
classical particles are of extremely small size. However, we may also see that an 
electron executing harmonic oscillation with frequency ujq has to accommodate 
relatively large values for x«, approximately equal to 0.01/y^Jo m - 

In summary, the probability axiom in the modified dBB quantum mechanics 
helps to distinguish the classical limit of quantum harmonic oscillator as one 
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Fig. 3. The complex trajectories in the n = 2 harmonic oscillator case. 
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Fig. 4. The complex trajectories in the n = 3 harmonic oscillator case. 

in which the oscillator is probable to be found only very close to the real axis. 
This result is very important for complex quantum trajectories, for it explains 
why the complex extension is not observable even indirectly in the classical 
limit. We anticipate that this property is true in other problems too. 



5 Conclusion 



The present work is a continuation of that in [T3] to obtain a probability den- 
sity for the entire complex x-plane, in the modified de Broglie-Bohm quantum 
mechanics. In the earlier work, a conserved probability (|6]) was proposed for 
particles in 1-dimensional stationary states, but only for those regions where 



trajectories enclose poles of the velocity field. This continues to be so in the 
present framework. But for those regions where trajectories do not enclose 
any poles, described as subnests, a conserved probability could not be found 
in [13] and it was considered a missing link in the formalism. We here adopt 
ty*(x)ty(x), which is equivalent to that in (J7|), as the extended probability den- 
sity for those special regions and note that the combined probability density 
helps to find answers to many pertinent questions that arise in this context. 

First of all, when compared to all other proposals for an extended probability 
density, the present combined probability has the important advantage that 
it is normalisable in the a^Xj-plane. At the same time, along the real axis in 
all regions, it agrees with the Born's probability, as required. But we note 
that the proposed probability for subnests does not obey a continuity equa- 
tion. However, it can be seen that for the entire subnested region, the net 
value of this probability is conserved. This property of local nonconservation 
of particles, also found in the FFM trajectory representation, is argued to be 
characteristic of quantum phenomena in general and is not to be dismissed as 
unphysical. In addition, it explains why the imaginary part of complex trajec- 
tories do not leave any observable imprints in the case of classical particles. 
The present scheme is such that for a classical harmonic oscillator, probability 
is considerable only for those trajectories which lie extremely close to the real 
line, whereas for an electron in harmonic motion, trajectories have substantial 
probability even when their imaginary part is appreciable. The observation 
that complex trajectories with large imaginary part are least probable in the 
classical limit is physically much relevant for the representation. 
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